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Q Abstract. It is widely believed that the heating of the chromosphere in 
quiet-Sun internetwork regions is provided by dissipation of acoustic waves that 
are excited by the convective motions close to the top of the convection zone 
fN) ' and in the photospheric overshoot layer. This view lately became challenged by 

observations suggesting that the acoustic energy flux into the chromosphere is 
too low, by a factor of at least ten. Based on a comparison of TRACE data with 
synthetic image sequences for a three-dimensional simulation extending from the 
top layers of the convection zone to the middle chromosphere, we come to the 
contradicting conclusion that the acoustic flux in the model provides sufficient 
energy for heating the solar chromosphere of internetwork regions. The role of 
a weak magnetic field and associated electric current sheets is also discussed. 
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1. Introduction 



In two recent papers iFossum h Carlssonl (2005, 2006, hereafter FC05 and FC06) 
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come to the conclusion that 'high-frequency acoustic waves are not sufficient to 
heat the solar chromosphere'. This conclusion is based on a study of image 
sequences from the 1600 A channel of TRACE (Transition Region and Coro- 
^ . nal Explorer). The continuum intensity at the wavelength of 160 nm originates 

from the upper photosphere. There, acoustic waves suffer no longer strong 
radiative damping and have not yet become strongly non-linear. Hence, the 
TRACE images seem to be an excellent means for measuring the total mechan- 
ical flux entering the chromosphere at its base. However, the mechanical flux 
cannot be measured directly but must be inferred from the measured inten- 
sity fluctuation (A///)i60nm- FC05 do this by comparing the average spectral 
power density for (A///)°gQ nm with those for corresponding synthetic intensity 

fluctuations (A///)^g u t ^ m from simul ations with the one-dim ensional radiation 
hydrodynamics code RADYN (e.g., ICarlsson Sz Steinl Il997l ). They find that 



the power spectra fit best for simulation runs that yield a mechanical flux at 
a height of 400 km (integrated over 5 to 50mHz) of merely 438 Wm -2 — ten 
times less than is needed to compensate the radiative losses of approximately 
4300 Wm " 2 (excluding L yman a) in the semi-empirical model of the chromo- 
sphere bv lVernazza et all (1981, hereafter VAL). 

Here, we present new ca lculations based on the three-dimensional radiation 
hydrodynamic simulation bv lWedemever et al.l (2004, hereafter W04), which ex- 
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hibit a dynamic pattern on a spatial scale too small to be resolved with TRACE. 
We discuss potential implications for the heating of the solar chromosphere. 



2. Acoustic heating on small spatial scales 

2.1. Three-dimensional non-magnetic simulations 

The three-dimensional radiation hydrodynam ics simulations by W 04 were car- 
ried out with the computer code C0 5 BOLD ( Frevtag et al.1 12002) and cover a 



horizontal region of 7 '.'7 x 7 '.'7 and extend from the upper convection zone at 
a depth of 1400 km to the middle chromosphere at a height of 1710 km. The 
spatial resolution is 40 km in the horizontal directions and 12 km vertically in 
the atmosphere. Magnetic fields are not included so that this model represents 
very quiet internetwork regions only. 

In order to keep the computations feasible, the radiative transfer is treated 
grey (i.e. frequency-independent) and under the assumption of local thermody- 
namic equilibrium (LTE). These are severe assumptions for the upper photo- 
sphere and the chromosphere above but are reasonable for the low photosphere 
in the framework of time-dependent three-dimensional simulations. The model 
chromosphere is - in sharp contrast to semi-empirical static models like, e.g., 
those by VAL - highly intermittent in time and space, consisting of a dynamic 
mesh-work like pattern of hot shock wave fronts and cool post-shock regions in 
between. The assumptions made for the radiative transfer affect the tempera- 
ture amplitudes and the chromospheric energy balance but are reasonable for the 
low photosphere, where the acoustic waves are excited. The model provides suf- 
ficient acoustic energy to counterbalance the radiative losses in the internetwork 
chromosphere w ithout need for fu rther contributions from processes related to 
magnetic fields ( Wedemeverl 120031 ) . 



2.2. Intensity synthesis 

The general procedure for our analysis can be summarised like this: (i) Synthesis 
of the continuum intensity at relevant ultra-violet (UV) wavelengths, (ii) integra- 
tion of the intensity over wavelength to a synthetic TRACE channel, (iii) image 
degradation to the spatial resolution due to the instrument's optics, (iv) spatial 
intensity integration for mimicking the TRACE detector (pixels), and finally 
(v) calculation of (average) power spectra from the intensity sequences resulting 
from step (i) to (iv). 

We use a short simulation sequence from W04 with a duration of lOmin 
and a cadence of 10 s. For each snapshot we calculate continuum intensity maps 
with the radiative transfer code Linfor3D (http://www.aip.de/~msteffen/) 
at wavelengths from 150 nm to 180 nm with an increment of 5nm, covering 
the wavelength range relevant for the TRACE 160 nm channel. With increas- 
ing wavelength the intensity fluctuations become smaller, whereas the absolute 
intensity level grows strongly. Including longer wavelengths thus effectively re- 
duces the signal I / < I > and the resulting spectral power density. The inten- 
sity sequences for different wavelengths (150 nm - 180 nm) are now in tegrated 
according to the instrument's transmission curve (iHandv et al.lll999l . see also 
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Figure 1. Synthetic intensity images for 160 nm passband of TRACE: 
a) original image, b) degraded to a spatial resolution of 0'/5, and d) af- 
ter integration on TRACE 0'.'5 pixels. The image degradation is repeated for 
a coarser resolution of l'.'O in panel c) but with the same pixel size in panel e). 
The dashed contours for the original intensity in panel a help to identify the 
enclosed regions in the other panels. Panel f) shows observed TRACE data 
from September 28 th , 2005, taken from a data set kindly provided by Fos- 
sum & Carlsson. The dashed squares mark pixels that would be disregarded 
according to the pixel selection method by FC06. The grey scale ranges of 
the panels are clipped individually. The maximum value is given in brackets 
above each legend. The intensity contrast is noted on the side of each panel. 
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http://www.lmsal.com/solarsoft), resulting in a short sequence for the syn- 
thetic TRACE channel. A typical intensity image is shown in Fig. QJi. 

In the next step the synthetic intensity images are convolved with a point- 
spread function (PSF) in order to mimic the instrument's UV optics (Fig. [lb). 
Unfortunately, the PSF of TRACE is poorly known but might be represented 
bes t by a Gaussian kernel of width (FWHM) l arger than 0^5 but not more than 
1" ([Kriiger et al.ll200ll : [Schneider et al.ll2004l : ISchriiven 120061 1. The detector of 
TRACE has a pixel-scale of 0.5 "/px. A synthetic pixel thus corresponds to 
spatial integration over (square) regions of 9 x 9 model grid cells (Fig. [lb). Note 
that FC05 additionally apply a binning of 2 x 2 pixels but FC06 use no binning. 
An example of observational data from FC06 is shown in Fig. [TU for comparison. 

Figure Q] clearly demonstrates that the pattern on small spatial scales is 
lost and only the brightest features are still visible in the synthetic image. Con- 
sequently the intensity contrast decreases from 81 % to 42 % for the displayed 
snapshot (89% and 33%, respectively, for the whole sequence). The lower in- 
tensity contrast of the observational data (35 % for the displayed region, 30 % 
for the whole time series for that region, and 26 % for the time series of all pixels 
according to the mask by FC06) suggests that the true resolution of TRACE 
might be coarser than the 0'.'5 shown in Fig. [lb, which represents a lower limit 
only. A resolution of 1", on the other hand, might indeed be considered an 
upper limit as it reduces the contrast of the synthetic images to 33 % for the 
displayed snapshot and only 26 % for the whole time series. 



2.3. Spectral power density 

We derive the temporal spectral power density (one-sided) for each spatial po- 
sition in the synthetic image sequences separately. The horizontally averaged 
power spectra are shown in Fig. [2] for the sequence with the original spatial 
resolution of the input model and also for the two sequences degraded to lower 
resolution (0'.'5 and 1 " ). For comparison we plot the power spectra that we cal- 
culate from the observed TRACE sequence from September 28 th , 2005 (FC05), 
(i) for the pixel sample of FC06 and (ii) for the small region shown in Fig. [TH. 
Image degradation obviously suppresses power on small spatial scales. The av- 
erage power is reduced by a factor of 10 at 12 mHz for the degraded sequence 
with a resolution of 0'/5, with respect to the original image sequence, whereas 
a resolution of 1" results in a factor of 27 for the same frequency (see Fig. [2]). 
The reduction in power increases with frequency up to ~ 20 mHz as higher fre- 
quencies are connected to smaller wavelengths (of acoustic waves) that are more 
severely affected by a limited resolution than large-scale fluctuations. The power 
when using a PSF of 1" width agrees well with the observational results in the 
range 3 mHz to ~ 12 mHz. 

The discrepancy at low frequencies is due to the short sequence duration, 
which results in a frequency resolution (see solid circles) that is too small to 
resolve the local minimum at v ~ 2.5 mHz. A longer synthetic image sequence 
at 160 nm with a duration of 76min and a cadence of 30 s, that we calculated 
in addition, shows results similar to those for the short sequence at 160 nm with 
respect to absolute magnitude of power and the relative effect of degrading the 
spatial resolution. The longer duration provides much finer frequency sampling 
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Figure 2. Average one-sided power density spectra for the short sequence 
of synthetic TRACE images (lines with dots) for original spatial resolution 
(solid), degradation to 0'.'5 (dashed) and 1" (dot-dashed) without pixel- 
binning. The solid dots mark the data points. The two other curves are 
calculated from observational data from September 28 th , 2005, in the same 
way as done for the synthetic sequences. The thick curve includes all pix- 
els according to the mask by FC06, whereas the thin one is for the region 
displayed in Fig. [TU only. 



so that the average power spectrum indeed exhibits a local minimum at low 
frequencies. 

The increasing difference at higher frequencies is most likely connected to 
the assumption of LTE that had to be made in order to keep the calculations 
computationally tractable. The assumption is validated by comparison with 
a more detailed calculation with the RH code (see, e.g., lUitenbroekl l2000l ) for 
eight 2D slices from the model for a wavelength of 160 nm, taking into account 
deviations from LTE. The rms fluctuations of the resulting intensity samples 
are 150.8% for Linfor3D (LTE), 199.3% for RH in LTE, and 134.9% for RH 
in NLTE. The difference for the LTE samples must be attributed to (i) the as- 
sumption of 2D for our application of RH in contrast to full 3D with Linfor3D, 
(ii) differences in opacity data (in particular concerning photoionisation cross- 
sections of silicon), and (iii) the neglect of continuum scattering in Linfor3D. The 
latter, however, is expected to be of minor importance as there are few contin- 
uum scattering processes that are relevant at wavelengths around 160 nm. More 
important is the correct treatment of bound-free transitions of silicon (taken 
into account in RH). They can be regarded as bound- free "two- level" resonance 
scatterers and effectively smooth out fluctuations. Nevertheless this qualitative 
comparison demonstrates that neglecting non-LTE effects tends to produce too 
large relative intensity fluctuations and with that an excess in average power for 
higher frequencies as visible in Fig. [2J 



In addition, a simulation sequence with non-grey radiative transfer (jWedemever 

2003) was used for intensity synthesis at 160 nm. The resulting power spectra at 



original and lowered resolution are very similar to the grey sequence used here. 
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Figure 3. Logarithmic current density, log | j | , in a vertical cross section (top 
panel) and in four horizontal cross sections in a depth of 1180 km below, and 
at heights of 90 km, 610 km, and 1310 km above the mean surface of optical 
depth unity from left to right, respectively. The arrows in the top panel 
indicate the magnetic field strength and direction. The dashed line indicates 
the position of the vertical section, [j] = 3 • 10 5 A/m 2 . 



3. Ohmic dissipation of chromospheric current sheets 
3.1. Three-dimensional MHD simulation 

We have also carried out a three-dimensional simulation including magnetic 
fields. The computational domain extends from a depth of 1400 km below the 
mean surface of optical depth unity to 1400 km above it well into the chromo- 
spheric layer s. The horizontal exte nt is only slightly smaller than for the model 
by W04. See ISchaffenberger et al.1 (2005, 2006) for details. 



The MHD simulation starts with a homogeneous, vertical, unipolar mag- 
netic field of a flux density of 10 G superposed on a previously computed, relaxed 
model of thermal convection. This flux density is thought to represent magneto- 
convection in a very quiet network-cell interior. The magnetic field is constrained 
to have vanishing horizontal components at the top and bottom boundary but 
lines of force can freely move in the horizontal direction. Because the top bound- 
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ary is located in the chromosphere, the magnetic field can freely expand with 
height through the photospheric layers into the more or less homogeneous chro- 
mospheric field, different from conventional simulations that extend to a height 
of typically 600 km only. 

A very common phenomenon in this simulation is the formation of a 'mag- 
netic canopy field' that extends in a more or less horizontal direction over ex- 
panding granules and between photospheric flux concentrations. The formation 
of such canopy fields proceeds by the action of the expanding flow above gran- 
ule centres. This flow transports 'shells' of horizontal magnetic field to the 
upper photosphere and lower chromosphere, where shells of different field direc- 
tions may be pushed close together, leading to a complicated network of current 
sheets in a height range from approximately 400 to 900 km. 

This network can be seen in Fig. [3] (top), which shows, for a typical snap- 
shot of the simulation, the logarithmic current density, log|j|, together with 
arrows indicating the magnetic field strength and direction. Figure [3] (bottom) 
shows from left to right log \j\ in four horizontal cross sections in a depth of 
1180 km below, and at heights of 90 km, 610 km, and 1310 km above the mean 
surface of optical depth unity. Higher up in the chromosphere (rightmost panel) , 
thin current sheets form along shock fronts, e.g., in the lower left corner near 
x = 1.4 Mm. 



3.2. Ohmic dissipation 

In this simulation we have not taken an explicit magnetic diffusion into account 
so that the effective electrical conductivity is determined by the inherent ar- 
tificial diffusion of the numerical scheme. Therefore, due to lack of a realistic 
electric conductivity, we here use molecular values when computing the ohmic 
dissipation of the chromospheric current sheets. Although the values might be 
orders of magnitude too high, they still can be employed for the following cal- 
culation, which gives a rough idea about the significance of ohmic dissipation 
for chromospheric heating. The typical current density in the height range from 
400 to 900km of Fig. [3] (top) is j = 0.03 Am -2 . The electrical conductivit y, cr. 



in th e photosphere and the lower chromosphere is about 10 to 100 A/ Vm (jStix 
20021 ) . Using this value, the ohmic dissipation is 



Pj = - j 2 « 0.03 2 « 10- 5 . . . 10- 4 W m~ 3 . (1) 

3 a 10 . . . 100 V ' 

When integrating over a height range of 500 km this heat deposition leads to an 
energy flux of 5 to 50 Wm~ 2 . This value is about two orders of magnitude short 
of being relevant for chromospheric heating. However, the amplitude and width 
of the current sheets in the simulation is determined by the spatial resolution of 
the computational grid rather than by the molecular conductivity. The effective 
conductivity in the simulation may easily be two orders of magnitude lower so 
that magnetic heating by ohmic dissipation must be seriously taken into account 
as a chromospheric heating agent. Further simulations, taking explicit ohmic 
diffusion into account will clarify this issue. 
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4. Discussion 



The photosphere of the model by W04 c an account for many aspects of obser - 
vations of internetwork regions (see, e.g., Leenaarts <fe Wedemever-Bohm 2005). 
Here we provide evidence that the UV intensity in the TRACE 160 nm channel 
is matched reasonably well, too. When degraded to the spatial resolution of the 
instrument, the synthetic intensity images presented here have similar contrasts 
and also result in average spectral power densities similar to the TRACE ob- 
servations. The power spectra exhibit an enhancement around periods of 5min 
and otherwise a decrease with frequency as it was found by FC05 and FC06 
both from observations and their ID simulations. 

At higher frequencies, however, we derive too much power from the synthetic 
image series. This discrepancy is most likely due to the LTE assumption made 
for the intensity synthesis with Linfor3D. The test calculations with the RH code 
show that taking into account deviations from LTE reduces the intensity contrast 
by ~ 60 % compared to the LTE case for a wavelength of 160 nm. The deviation 
from LTE at the shorter wavelengths is mostly due to the strong over-ionisation 
of Si I, which affects the continuum forming bound-free transitions. Similar to 
true continuum scattering, this effect smoothes out fluctuations, so that the 
LTE assumption used here actually tends to produce too large fluctuations and 
thus an excess in power. We argue that the derived reduction of ~ 60 % in 
intensity contrast can be considered as an upper limit for the influence of this 
scattering-like process. In contrast, the deviations from LTE are much smaller 
for the longer wavelengths. For a rough estimate we derive an average ratio 
of LTE and NLTE intensities based on the RH calculations and use it to scale 
the LTE intensities of the whole synthetic 160 nm sequence. This "pseudo- 
NLTE" time series leads to a further reduction in power by approximately a 
factor of two around periods close to 5min to four for higher frequencies. The 
tendency of larger power reduction with increasing frequency might be even 
more enhanced by non-LTE effects being more important on time scales too 
short for relaxation to equilibrium conditions. The "pseudo-NLTE" estimate is 
of course an order of magnitude approximation only and needs to be replaced 
with detailed NLTE calculations that include not only 160 nm but the whole 
wavelength range relevant for the TRACE channel. Taking into account line 
blanketing acts in the opposite direction by reducing the intensity at the longer 
wavelengths and thus increasing the relative fluctuations. A comparison of the 
spectra calculated with MULTI for the VAL model A results in a reduction of 
integrated intensity by a factor of ~ 2 ( Carlsson 120061 ). But whether it really 
counterbalances the expected NLTE effects and if so by how much and for which 
frequencies remains subject to detailed radiative transfer calculations. 

The computational box of the hydrodynamical simulation is small (5.6 Mm 
x 5.6 Mm). Consequently, the number of possible oscillation modes is much 
smaller than in the real Sun so that the power in the model is distributed over 
a smaller number of modes. The total energy content of the modes, however, 
remains the same. Also the finite grid resolution limits the treatment of high- 
frequency waves but FC06 state that the neglect of frequency beyond 20 mHz is 
not an important omission. 

When computing the mechanical energy flux from the intensity fluctuations, 
FC05 and FC06 assume that the resolution of TRACE is sufficient to fully cap- 
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ture the energy carrying waves. FC0 6 state th at tests with 2D models related 
to the hydrodynamic simulations by iHansteenl (2004) imply that only a factor 
of two in acoustic flux is missed due to insufficient spatial resolution. In con- 
trast our 3D study presented here shows that the unresolved small-scale pattern 
provides more power by an order of magnitude or even more. 

At this stage we cannot exactly determine how much of the additional 
power hidden on small spatial scales is indeed transported to the chromosphere 
and how much remains in the photosphere and therefore does not contribute to 
chromospheric heating. Separating the acoustic flux from the total mechanical 
flux is not easy for our 3D model the latter gets on average even negative in 
the middle photosphere as a consequence of convective overshooting. There 
is a local maximum of the radiative flux divergence in the low chromosphere 
and a mechanical counterpart of opposite sign that we interpret as evidence for 
the conversion of mechanical energy provided by acoustic waves into therma l 
energy, which is subsequently emitted in form of radiation ( Wedemev er 2003). 
Integration of the mechanical flux divergence over the height range from 570 km 
to 1710 km results in an acoustic energy flux of 7.9 kWm" 2 which is sufficient to 
counterbalance the (semi-)empirically determined radiative energy losses from 
the chromosphere (see, e.g, VAL). Moreover, the total mechanical flux at a 
height of 800 km (which approaches the pure accousti c flux there) amounts to 
3.8 kWm~ 2 and is thus in line with the ID models by iRammacherl (2005). 

On the other hand, the energy flux due to Ohmic dissipation of current 
sheets in the MHD model chromosphere is only of the order of 5 to 50 Wm~ 2 . 
This value must be considered as a lower limit only because of the uncertainty 
in the electrical conductivity. It can thus not be ruled out that the magnetic 
heating contribution might be underestimated by even two orders of magnitude. 

Although the energy balance of the model chromosphere is affected by 
a too simple treatment of the radiative transfer, the calculations presented 
here suggest that the acoustic contribution to chromospheric heating cannot 
be neglected. Magnetic fields may a lso play an indirect role for guiding acous 
tic waves and for mode conversion (jSuematsul[i990l : lUlmschneider et al.l fl99ll 
Jefferies et aDl2006l 1. 



5. Conclusions 

Based on the analysis presented here we conclude that TRACE might miss 
acoustic power by at least one order of magnitude due to its limited spatial reso- 
lution. Horizontal flows, as they are clearly seen within the reversed granulation 
pattern, contribute to the intensity variations although their relevance for the 
energy transport to the upper layers is not clear yet. 

Matching the empirical power spectra with a model that has enough acous- 
tic power to counterbalance the observed chromospheric emission should not 
be interpreted such that high-frequency waves alone are sufficient to heat the 
chromosphere. It is not clear from this study how important the small rela- 
tive power contribution of high-frequency waves is for chromospheric heating. 
Moreover, the model used here exhibits a chromosphere with pronounced dy- 
namics and a highly inhomogeneous structure - in sharp contrast to static and 
one-dimensional semi-empirical VAL-like models. The ubiquitous shock waves, 
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which produce our model chromosphere's characteristic structure, play certainly 
an important role for the energy balance of our model chromosphere. 

It must be emphasised that - owing to simplifications that are necessary 
in order to keep the problem computationally tractable - the topic should be 
addressed again with detailed NLTE calculations. Nevertheless, the presented 
study clearly demonstrates that three-dimensional radiation hydrodynamics sim- 
ulations in principle can reproduce the TRACE observations of internetwork 
regions without strong need for magnetic fields, although weak fields are most 
likely an integral par t of those regions. We conclude that an acoustic "basal 
flux" ( Schrijver 19871 ) must still be considered a possible option for heating the 



chromosphere of internetwork regions. 
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